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ABSTRACT 


Massive stars can shed material via steady, line-driven winds, eruptive outflows, or mass-transfer 
onto a binary companion. In the case of single stars, the mass is deposited by the stellar wind into the 
nearby environment. After the massive star explodes, the stellar ejecta interact with this circumstellar 
material (CSM), often-times resulting in bright X-ray line emission from both the shock-heated CSM 
and ejecta. The amount of material lost by the progenitor, the mass of ejecta, and its energetics all 
impact the bulk spectral characteristics of this X-ray emission. Here we present a grid of core-collapse 
supernova remnant models derived from models for massive stars with zero age main sequence masses 
of ~ 10 — 30 Mo evolved from the pre-main sequence stage with wind-driven mass-loss. Evolution is 
handled by a multi-stage pipeline of software packages. First, we use mesa (Modules for Experiments in 
Stellar Astrophysics) to evolve the progenitors from pre-main sequence to iron core collapse. We then 
use the Supernova Explosion Code (snec) to explode the mesa models, and follow them for the first 100 
days following core-collapse. Finally, we couple the snec output, along with the CSM generated from 
mesa mass-loss rates, into the Cosmic-Ray Hydrodynamics code (ChN) to model the remnant phase 
to 7000 years post core-collapse. At the end of each stage, we compare our outputs with those found 
in the literature, and we examine any qualitative and quantitative differences in the bulk properties of 
the remnants and their spectra based on the initial progenitor mass, as well as mass-loss history. 
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1. INTRODUCTION 


Core-collapse (CC) supernova remnants (SNR) are 
important drivers of galactic evolution, as well as tracers 
of massive star populations and their evolution. Broad- 
band observations have yielded a wealth of information 
about remnants, but coupling observations to theory 
requires sifting through a great deal of uncertainty in 
the parameters that dictate evolution (e.g., Patnaude 
& Badenes 2017). In particular, stellar mass loss rates 
throughout the life of the progenitor star can have a 
strong impact on the structure and dynamics of the re- 
sulting remnant (Patnaude et al. 2015, 2017). 

Many open questions remain concerning the efficiency 
of both steady and episodic mass loss mechanisms, the 
physical conditions in the progenitor required to drive 
mass loss, and exactly when in the stellar life cycle 
certain mass loss mechanisms are more prevalent (e.g., 
Quataert & Shiode 2012; Smith 2014; Fuller 2017). In 
the case of some Type IIP SNe, radio observations in- 
dicate quasi-steady mass loss up to nearly the point of 
core collapse (e.g., SN 2011ja; Chakraborti et al. 2016), 
while optical observations of so-called "supernova im- 
postors" suggest that some progenitors undergo mas- 
sive, eruptive events only a few years before the onset of 
core collapse (e.g., SN 2009fk; SN 2009ip; Moriya et al. 
2011; Margutti et al. 2014). Thus, current theories of 
stellar mass loss focus on wind-driven versus episodic 
events following the onset of core-carbon burning (Pat- 
naude et al. 2017). An important probe of these mass 
loss mechanisms comes from the forward shock emission 
of young SNR. The forward shock at early post-core- 
collapse time interacts with the CSM sculpted by late- 
time mass loss (Chevalier 1982a). Forward shock evolu- 
tion can therefore probe CSM density and composition 
at increasing radii around the SNR, and these results 
can yield the mass loss rates and durations as a func- 
tion of age pre-core collapse (Chevalier 1982b; Chevalier 
& Fransson 1994). 

While mass loss rates can depend strongly upon the 
mass loss mechanism, they also depend on parame- 
ters of the progenitor star such as the zero-age main- 
sequence (ZAMS) mass, the stellar metallicity, rotation, 
and whether the star is in a binary system (e.g. Maeda 
et al. 2015; Ouchi & Maeda 2017). Overall, the complex- 
ity of simply incorporating all these properties into a 
realistic model for the progenitor is an intractable prob- 
lem. Likewise, fully realized multidimensional simula- 
tions of core collapse supernovae through and beyond 


shock breakout are not readily available and are often 
tailored to specific objects (e.g., Orlando et al. 2020, 
2021), nor are evolutionary models for supernova rem- 
nants in complex circumstellar environments. Finally, 
small refinements to model input parameters derived 
from observations can render a multidimensional sim- 
ulation obsolete, and any one high-fidelity end to end 
model may not be applicable to the diverse types of core 
collapse supernovae and supernova remnants we observe 
in the local Universe. 

To address the complexities introduced when trying 
to develop a high fidelity model for the end-to-end evo- 
lution of a massive star from the pre-main sequence 
through the SNR phase, we present here a grid of one- 
dimensional (1D) models for core collapse SNR evolu- 
tion, derived from models for massive stars. We use 
the "Modules for Experiments in Stellar Astrophysics" 
(mesa; Paxton et al. 2011, 2013, 2015, 2018, 2019) for 
the stellar evolution, assuming purely line driven winds 
with rates from van Loon et al. (2005) and de Jager et al. 
(1988). At the onset of core collapse, the progenitor 
models are coupled to the "Supernova Explosion Code" 
(snec; Morozova et al. 2015) modified to include ex- 
plosive nucleosynthesis (Patnaude et al. 2017). Finally, 
the output from snec is coupled to our "cosmic-ray hy- 
drodynamics" code (ChN; Ellison et al. 2007; Patnaude 
et al. 2009, 2010; Lee et al. 2014), to simulate the SNR 
evolution. 

We focus on producing a wide range of progenitors 
such that we can cover a reasonable amount of stel- 
lar mass parameter space in 82.1, and calibrate our SN 
models such that they are in reasonable agreement with 
simulated compositions from the literature in 82.2. Fi- 
nally, we examine the remnants as a function of time 
after core-collapse with particular attention paid to X- 
ray spectral lines for higher Z elements (eg. S and Fe) 
in 82.3. We examine centroid information to determine 
the behavior of observables for the current generation 
of telescopes in 83.1. We then move to a discussion of 
He-like emission from high Z elements and discuss what 
additional power these probes may have in determining 
progenitor information in 83.2. Additionally, we apply 
our models to the well-studied Galactic SNRs in 8 4. 
We close with a discussion of what spectral information 
will be available to future X-ray missions, and discuss 
the behavior of some observables for spatially unresolved 
remnants in 85. 
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2. METHODS 


Our current grid of models encompasses ~ 570 mod- 
els with ZAMS masses ranging from 9.6Mo to 30.0Mo. 
Each progenitor mass was modeled using various pa- 
rameter settings at each stage of the pipeline, and the 
results are contained within publicly available git repos- 
itories for ease of versioning and to aid in reproducing 
our results. In this section, we undertake a detailed de- 
scription of the pipeline, as well as the various parameter 
choices used to generate the final models. 


2.1. Progenitor Modeling in mesa 


We begin our pipeline with the stellar evolution code 
mesa using version mesa-r10398 (Paxton et al. 2011, 
2013, 2015, 2018, 2019). Models were derived from the 
sample project example make pre ccsne. We chose to 
use a pre-existing set of inlists for a number of reasons, 
the chief ones being ease of implementation and stabil- 
ity across a wide range of physical parameters, as well as 
ease of reconstructing our pipeline against future mesa 
versions. For our baseline models, we initially left all 
parameters at their default, aside from Z which we in- 
creased to Zo. This resulted in non-rotating stars con- 
figured to use the 'Dutch' wind scheme for mass loss in 
the red supergiant (RSG) phase, with rates derived from 
(de Jager et al. 1988). These settings were applied to a 
grid with ZAMS masses ranging from 9.0Mo X Mzams 
< 30M. Of the 205 models, ~ 90% successfully reached 
core-collapse, with the remainder encountering difficulty 
during the Si burning phase. These failures were biased 
toward the lower end of the mass range (~ 9M), with 
no stars below 9.6Mọ reaching core-collapse. This is 
unsurprising as neutron-rich material plays a larger role 
in the later burning stages, and the 21-isotope network 
used does not contain enough of these isotopes to accu- 
rately model these lower-mass stars (Paxton et al. 2015). 

Following our non-rotating models, we performed an 
additional run over the same mass range with a speci- 
fied angular momentum of 2 = 0.550 crit, and an initial 
surface velocity of v = 0.55verit- Here verit and OQerit are 
defined as the maximum values of the surface velocity 
and total angular momentum that can be maintained 
by the model star before it becomes gravitationally un- 
bound. These initial conditions were chosen to match 
with those used by Renzo et al. (2017) in their explod- 
ability study. The final masses of the models including 
the compact object are plotted against the ZAMS mass 
in Figure 1 for the rotating de Jager models. In spite 
of the large initial mass range, the final mass range en- 
compasses a narrow set of masses, with the average final 
mass being ~ 12.5Mo. 
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Figure 1. Final Mass vs. ZAMS mass for the rotating de 
Jager wind models. Shaded region corresponds to Mena = 
12.3977 65 Mo. Although ZAMS masses span ~ 20 Mo, the 
majority of final masses are confined to < 14Mo. 


In examining this initial set of models, we discovered 
several that contained a large loop in their HR diagrams 
of constant luminosity, but increased effective surface 
temperature (Teg), which ultimately returned to the ex- 
pected late-stage evolution seen in the non-looping mod- 
els. These loops appeared to be blue-loops, as seen in 
the evolution of many stars (e.g. Xu & Li 2004; Paxton 
et al. 2019). Unlike observed blue loops, the modeled 
stars experienced significantly reduced mass-loss due to 
the modeled hot winds being less effective at lifting mass 
from the star than the cool winds of stars that did not 
exhibit a loop. A closer examination showed that these 
stars experienced the same Fe opacity peak that would 
be expected for looping stars, but it did not explain why 
others did not experience this (Paxton et al. 2019). Per- 
forming a mesh resolution study seemed to indicate that 
these loops are physical, as although the stars that orig- 
inally exhibited the loops may have lost them at higher 
resolutions, the looping behavior was a constant, and 
ultimately settled to a narrow selection of masses. 

In addition to these models, we have two secondary 
grids containing models from 12 — 23Mọ. One was run 
with a cool wind based on the van Loon wind scheme 
(van Loon et al. 2005), and the other with the metallic- 
ity set to the default Z of 0.006. Combined, these groups 
give us a total of 570 models simulated from pre-main- 
sequence to iron core-collapse. The final models were 
recorded along with data for several key outputs during 
the lifetime of the star, the chief ones being mass loss, 
wind velocity, temperature and luminosity. The former 
two parameters are important for determining the be- 
havior of the circumstellar environment, the latter two 
are important for determining where in the stellar life- 
time certain mass loss events or significant changes have 
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Figure 2. Example HR diagrams for four progenitor models, 
exhibiting the three major evolutionary tracks. We note that 
three of the four models explode as RSG stars, while the most 
luminous star at CC explodes as BSG. 


occurred, as well as for making sure the models are well 
behaved throughout the process. We present example 
HR diagrams for the major evolutionary tracks in Fig- 
ure 2, and discuss any variation based on the parame- 
ters. 


O 


> 


Pinal Effective Temperature (K) 


5.0 17.5 20.0 22.5 25.0 27.5 30.0 
Initial Mass (Mo) 

Figure 3. Final Temperature vs. ZAMS mass for the rotat- 
ing models (black), and the van Loon wind models (grey). 
The dashed line corresponds to the hottest star in the van 
Loon interpolation range. The dotted line indicates solar 
Teg. Stars with a ZAMS mass below ~ 23 Mo result in RSG 
progenitors (hatched region). Above that point, stars that 
do not experience drops in mass-loss in the carbon-burning 
stages begin to explode as hotter and more compact objects. 


The cutoff for RSG progenitors occurs ~ 22 Mo. 
Beyond that point, the effective temperature increases 
rapidly for more massive progenitors, and they quickly 
become yellow supergiant (YSG) and blue supergiant 
(BSG) stars, as can be seen in Figure 3. We note this 
is broadly consistent with Katsuda et al. (2018), who 


found that observed RSG progenitors are similarly lim- 
ited. This is further born out by Davies & Beasor (2018) 
who found that when considering non-RSG progenitors, 
along with other types of CCSNe, there is no strong 
indication of a mass limit on progenitors in general. 
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Figure 4. Mass loss rate parameter space for the rotating de 
Jager models. The parameter space is bounded by a 30M 
star above, and a 10M star below. Note how the overall 
trend exhibits a significant ramping up of mass loss in the last 
500 — 800 kyr. Higher mass stars exhibit an earlier increase 
around ~ 1 — 1.3 Myr before CC, but the 800 kyr peak 
removes nearly an order of magnitude more mass. The low 
mass stars exhibit very low mass-loss rates across the board, 
so it becomes less important to determine the exact point 
where mass-loss reaches RSG levels. 


We can similarly dissect the mass-loss history of the 
star to determine what points would be most important 
for modeling the CSM that would be close enough to 
interact with the subsequent remnant. In theory, we 
can model the entire wind from ZAMS to core-collapse, 
but in practice, the main-sequence winds have already 
reached = 15pc by core-collapse, and would thus be un- 
likely to interact with the remnant during the simulated 
timeframe. As such, we focus on the material removed 
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from the star during the RSG phase, which generally 
occurs during the final 10° — 10° years. This becomes 
readily apparent when examining the difference in mass- 
loss rates during the main sequence and early post-main 
sequence versus the final years of stellar life. The range 
of the mass-loss histories over the last two million years 
of evolution is presented in Figure 4. 

As mentioned, many of our models reach CC as RSG 
stars. With this in mind, we examined the impact a 
dust-driven wind as described by van Loon et al. (2005) 
would have on mass-loss rates. The change in mass-loss 
can be significant, especially for lower mass stars, as can 
be seen in Figure 5, but the applicability of the van Loon 
wind scheme is limited for our models. Even our coolest 
stars are at the upper end of the interpolation range for 
the wind prescription (Figure 2), and these still spend a 
large fraction of their post-main-sequence life above the 
required temperature. Even with these caveats though, 
van Loon RSG winds are in reasonable agreement with 
the derived mass-loss rates (M) for 2002hh, 1999em, 
and 2004et, which all have derived masses within the 
van Loon interpolation range of our models. At higher 
masses, the van Loon wind does not cover much of the 
post main-sequence lifetime, and the increase in mass 
loss is only slightly more than expected in the rotating 
model case. 
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Figure 5. The average mass loss rates for the final 500 kyrs 
before CC as a function of ZAMS mass for the rotating de 
Jager models (blue line), the non-rotating de Jager models 
(orange line), and the non-rotating van Loon models (grey 
line). We compare these to derived mass-loss rates and their 
uncertainty regions (colored boxes) for IIp SNe as presented 
in Dwarkadas (2014). Many SNe align with the model mass- 
loss rates. 


As shown in Figure 5, all our models have wind- 
driven mass-loss rates of M ~ 1077 — 107? M /yr, which 
is consistent with previous studies (e.g., Dwarkadas 
2014). There is evidence of larger mass-loss rates (2 


^2 


10-4Mo/yr) for type II SNe in the final years before 
CC (e.g. Forster et al. 2018). The exact affect this has 
on the CSM will be discussed in 2.3, but it would not be 
noticeable in our averaged mass-loss rates due to their 
short lifetime compared to the total post-main sequence 
period. The total mass loss was dictated by the chosen 
wind mechanism with the parameters set identically for 
all models. For stars at the lower end of the mass range, 
the RSG winds are not sufficient to remove a significant 
fraction of the hydrogen envelope, as seen in Figure 6. 
The H envelope mass is also well correlated with the 
final effective temperature. Based on the H-envelope 
and mass-loss data, partially stripped envelope SNe (ie. 
IIL) do appear to be possible above M ~ 20Mo, but 
lower masses will need additional mass-loss mechanisms 
to produce these types of SNe. Stars above ~ 25M can 
produce IIb SNe based on Nomoto et al. (1995), but Ib/c 
SNe will also need additional mass-loss mechanisms at 
all masses (e.g., Sravan et al. 2019). 
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Figure 6. The Final Hydrogen Envelope mass as a function 
of ZAMS mass. Much like in Figure 3, The H-envelope size 
is more or less constant until M ~ 20Mo, at which point the 
envelope mass shrinks significantly. SNe type designations 
are based on Nomoto et al. (1995). The light grey points 
correspond to stars that have a Teg above the interpolation 
range for the van Loon wind scheme. 


The carbon-oxygen (CO) core is strongly tied to the 
ZAMS mass, and seems largely insensitive to the mass- 
loss rates of individual models, or to the mass-loss pre- 
scription. This is particularly obvious in Figure 7 as 
the outliers seen in Figures 5 and 6 are absent. This 
points to the CO core mass being a useful proxy for 
determining the ZAMS mass of a remnant as stated in 
Katsuda et al. (2018). Figure 8 presents three example 
models with similar final masses. The growth of the CO 
core and the loss of the H-envelope are readily apparent. 
'The abundance peaks are due to lighter elements burn- 
ing and mixing inward at the shell boundaries. Because 
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Figure 7. The Final Carbon-Oxygen Core mass as a func- 
tion of ZAMS mass. The CO-core is insensitive to mass- 
loss, unlike many of the other metrics described. We have 
included example data from Katsuda et al. (2018), which in- 
dicate the major mass bins they determined for the CO and 
ZAMS masses, using the measured Fe/Si mass ratio. 


the final compositions are not strongly affected by the 
mass-loss rates, and the final mass range is quite narrow, 
we choose to focus our discussion on the rotating rem- 
nant models and note that non-rotating van Loon (de 
Jager) models largely lead to an increased (decreased) 
CSM density, and otherwise follow the general trends 
discussed below. 


2.2. Modeling Supernovae with snec 


When the mesa models reach core-collapse, deter- 
mined to be when Fe core infall exceeds 1000 km s~!, we 
save the final model and pass it to the Supernova Ex- 
plosion Code (snec). As discussed in Patnaude et al. 
(2017), snec (Morozova et al. 2015) has been modi- 
fied to include explosive nucleosynthesis (the so-called 
approx21 network), and a Helmholtz Equation of State 
(Timmes & Swesty 2000). The progenitor models are 
exploded with snec, and evolved to an age of 100 days. 
The end product of the snec modeling is a model for 
the dynamics and composition of the ejecta, which can 
be read into ChN (see § 2.3). 

Besides the input progenitor model, snec allows the 
user to vary the explosion energy, the energy injection 
rate, the mass over which the energy is spread, and the 
neutron star mass cut. By tuning these parameters, the 
final chemical yields can be changed in the regions of 
the model undergoing explosive burning. Because of this 
unbounded parameter space, we chose to calibrate our 
models against published yields from one-dimensional 
SNe models. In particular, we benchmarked our models 
against Young & Fryer (2007, hereafter YF07), though 
a similar calibration against other abundance sets could 
be explored. 


Fractional Abundance (Xj) 


4 6 8 10 12 
Mass Coordinate (Mo) 


Figure 8. Stellar abundances as a function of Mass 
coordinate for the final mesa model for Mzams = 
(17.5, 20.6, 28.9) Mo (bottom to top). These three masses 
were picked as they all have a final mass ~ M. The growth 
of the CO core, and the loss of the H-envelope are readily ap- 
parent. Additionally, there emerges an independent Si core 
in the 29Mọ model. The peak Si abundance occurs within 
the CO core in the lower mass models. The vertical black 
line indicates the location of the excised mass for the central 
compact object. 


2.2.1. snec Explosion Calibration 


YFO7 use a 23M, ZAMS mass star as their baseline 
model, and they also implement a much larger burning 
network than our models which allows for many more 
element creation and destruction pathways. The prac- 
tical impact of using a larger network is in how much 
of each element is burned in a particular T — p regime. 
For instance, in the approx21 network, the photodisin- 
tegration reaction ??Ne(», o)!9O reaction (Arnett 1974) 
enhances 74Mg by a factor of X 10 over networks with 
larger numbers of isotopes. The reason for this excess 
is that in the more complex networks, neutrino cooling 
in the burning layers leads to increased densities which 
can prevent vigorous Ne ignition (Farmer et al. 2016). 

Similar enhancements exist for ?5Si burning, which 
burns to form Fe-group elements, including ?9Ni. The 
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Figure 9. Supernova Elemental Abundances plotted against 
Results from Young & Fryer (2007). The blue lines corre- 
spond to the yields for the SNe grid. Our abundances differ 
from YFO07 in part because of our much smaller reaction net- 
work, and lack of reaction chains containing neutron rich 
elements. 


approx21 network includes a "neutron-eating" isotope 
of chromium. Owing to a large Coulomb barrier, ?5Si 
does not generally burn directly to °°Ni. Instead, the 
reaction proceeds as first ?°Si photo-disintegration and 
then a chain of reactions in quasistatical equilibrium 
(QSE), resulting in Fe-group elements (Hix & Thiele- 
mann 1996). QSE reactions lead to a high degree of 
neutron-rich isotopes. However in the reduced reaction 
network we use here, the smaller number of isotopes 
used leads to a smoother composition profile, as seen in 
Figure 8. In practical terms, this can impact the syn- 
thesized X-ray spectrum of the final remnant. 

The chief difficulty in calibration was to match the Ni 
yields following the end of nuclear burning in the SNe. 
We used the explosion parameters in YFO7 as a starting 
point, and varied the injection time and injection region 
as well as the delay to align our final Ni yield with that 
in the literature. We present the yields of the surface 
rotating models plotted against the yields from YFO7 in 
Figure 9, and present the 23M calibration models in 
Figure 10. Most elements align reasonably well with 
our models, although the abundances for some mod- 
els are significantly higher than expected. ?4Mg over- 
shoots in part because the reaction network, approx21 
produces more Mg through longer Ne burning (Farmer 
et al. 2016). Additionally, as mentioned above, some 
isotopes such as !^N are burned in a limited number of 
reactions compared to larger nets (Timmes 1999; Paxton 
et al. 2011). 


2.3. Modeling Supernova Remnants with CSM 
interactions using ChN 
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Figure 10. Supernova Elemental Abundances plotted 
against Results from Young & Fryer (2007) for a 23Mo ex- 
ploded using various snec input paramters. Varying the ex- 
plosion parameters had a more pronounced effect on higher-Z 
elements. 


Once the supernovae have been simulated out beyond 
the point of Ni+Co decay (~ 100 days), the result- 
ing profile is piped into the Cosmic-Ray Hydrodynamics 
Code (ChN). ChN is a one dimensional Lagrangian hydro- 
dynamics code developed out of VH1 (Colella & Wood- 
ward 1984). The code simulates the dynamics of the 
remnant with initial conditions from the final composi- 
tion and energies supplied by snec, and self-consistently 
determines the ionization fractions of shocked mate- 
rial, as well as the thermal and non-thermal spectrum, 
using self-consistent non-equilibrium ionization calcula- 
tions (Patnaude et al. 2009, 2010; Lee et al. 2014). ChN 
output is coupled with atomdb and used to generate the 
broadband line emission from the remnant (Smith et al. 
2001; Foster et al. 2017). ChN additionally has the ability 
to simulate the effects of cosmic-ray acceleration on the 
remnant dynamics, but we defer an conversation those 
by operating in the so-called test particle limit (Ellison 
et al. 2007). We performed simulations to tsnr = 7000 
yrs. Similar work has already been performed for type 
Ia remnants, and the results have shown ambient den- 
sity driven differences in the evolution of the remnants 
(Martínez-Rodríguez et al. 2018). 

Initial plans also included modeling the dynamics of 
the CSM using VH1 and using that final profile as the ini- 
tial CSM in ChN. In practice, these simulations resulted 
in main-sequence bubbles several parsecs from the star 
and what was left within range of the SN shock was 
not significantly different from taking the average mass- 
loss rate over the last 500 kiloyears before core-collapse. 
This is not necessarily true of the stars that exhibited 
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Figure 11. Example Forward Shock (FS) spectra for 
Mzams = (17.5, 20.6, 28.9) Mo at 1000 yrs post-CC (bot- 
tom to top). Higher energy emission is particularly depen- 
dent on the wind density, with the He-like Fe (dotted) peaks 
growing significantly as the density increases. The complex 
between 6.4 keV and He-like Fe is composed of emission from 
lower Fe ionization states (~ Fe XVII upwards) and their 
satellite lines. H-like Fe (Ly a, dashed) also grows and de- 
velops satellite lines with increased density. 


blue loops and had extended periods of lower mass loss, 
but the effects caused by those changes are not cur- 
rently considered. Wind velocities from mesa seemed 
to be systematically under-reported due to issues defin- 
ing the exact boundary of the photosphere. To combat 
this, we instead performed two runs for each model with 
Uwind = (5,15) km s^!. Example spectra for the three 
models from Figure 8 are presented in Figure 11. 


3. ANALYSIS 


In spite of the fact that the models shown in Figures 8 
and 11 have approximately the same final mass, their 
spectra are noticeably different. One reason for this is 
that most SNe and SNR are found to evolve into a range 
of density profiles which are dictated by the progenitor 


mass-loss rate and the wind velocity, and the X-ray emis- 
sion is dependent upon the circumstellar density. This is 
shown in Figure 12 where we use the inferred mass-loss 
rates from our models to construct the circumstellar pro- 
file around the progenitor at the time of core-collapse, 
assuming a wind velocity of 10 kms^!. We also include 
inferred mass loss rates for a selection of SN IIP/IIL and 
the Galactic SNR G292.04-1.8. 
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Radius (cm) 
Figure 12. CSM Density profiles for the SNR at 


Core-Collapse. The wind velocity is taken to be 
Vwind = 10 kms-^!, making the density dependent on only 
the model mass-loss rate. Derived CSM densities for various 
remnants assuming the same RSG wind velocity. The data 
are in reasonably good agreement with the model densities. 
SNe data are taken from Chakraborti et al. (2013, 2016) for 
SN 2011ja and SN 2013ej, Black et al. (2017) for SN 2013by, 
and Yadav et al. (2014) for SN 2012aw. The inferred mass 
loss rate for G292.0+1.8 is taken from Lee et al. (2010). 


As seen in Figure 12, increasing progenitor mass, in 
broad terms, maps to increased mass loss rates dur- 
ing the later stages of evolution. Mass-loss rates in 
SN IIP/IIL and some Galactic SNRs are inferred from 
radio and X-ray observations, and are in broad agree- 
ment with the rates used here. For the larger rates dis- 
cussed in Fórster et al. (2018), The enhanced mass loss 
would lie within 10!^cm, and would have a larger affect 
during the very early stages of SN evolution, whereas 
the RSG winds drive the SN and SNR evolution at later 
epochs. Additionally, the red supergiant wind will inter- 
act with the low density wind from the main sequence 
star. The RSG wind sweeps the main sequence wind 
material into a radiatively cooled shell located ~ 107° 
cm from the progenitor (Castor et al. 1975). Models 
which extend into the boundary region between the red 
supergiant and main sequence winds will also need to 
consider radiative losses, as the swept up shell of ma- 
terial can rapidly decelerate the blastwave. This is an 
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important consideration for remnants with ages > 104 
yrs, but we do not consider it here. 


3.1. Calculating Centroid Energies and Ionization 
States 


X-ray spectra in supernova remnants encode physics 
about the chemical composition of the ejecta and the 
evolution of the progenitor. Line emission is produced 
by ionization and excitation of shocked circumstellar 
and ejecta material. In the case of SNRs, Ionization 
and recombination rates are driven by Coulomb inter- 
action, which are directly proportional to the CSM den- 
sity. Photoionization is not currently considered. The 
ionization state of each element is also driven by electron 
number, with higher-Z elements requiring more energy 
to fully ionize, while more readily recombining as can be 
seen in Figure 13. Therefore, line emission from these 
highly-ionized high Z elements strongly tracks the evo- 
lution of the remnant on timescales suitable for probing 
the late phases of stellar evolution. 

In order to better examine the trends that result from 
the CSM, Figure 13 and later figures that connect ion- 
ization to progenitor mass are plotted against the mass 
ratio 


M — Mas 
Mass Ratio = -74 model excised (1) 


Mzams 


Unlike the ZAMS mass, the mass ratio is sensitive to the 
mass-loss rate of the individual stars, so stars with rel- 
atively low mass-loss rates (ie. the blue loops) and low 
CSM densities will be grouped together. Looking at the 
ratio also allows for the larger final progenitor mass to 
be included in these analyses as the mass-loss rates for 
many blue loop stars are comparable to lower mass ob- 
jects but the ratio will be higher due to due to the much 
larger ejecta mass in the blue loop case. The differences 
in remnant composition do have implications for the ion- 
ization and spectral behavior, so it is important to look 
at metrics that consider all of these features. 

Higher ionization states of Fe have proven to be use- 
ful probes of remnant evolution due to the relatively 
high elemental abundance in the remnant and the rela- 
tive isolation of the K-shell emission in the X-ray spec- 
trum. Fe-Ka centroids have also been cited as a prob- 
able discriminant of SNR progenitor type (Yamaguchi 
et al. 2014). Core-collapse remnants show a strong trend 
toward the He-like transitions, while Type Ia centroids 
are dominated by Fe XVII and other similar ionization 
states. Our models support the general trend of core- 
collapse remnants exhibiting harder Fe-Ka centroid en- 
ergies, with models trending toward the Fe XXV reso- 
nance transition as the remnant ages. We do note that 
the circumstellar density is the strongest determinant of 
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Figure 13. Average ionization states in the shocked re- 
gion for the elements Fe, S, and Si (top to bottom) in a 
fast wind (solid, vina = 15 kms~'), and a slow wind 
(dashed, Vwing = 5 kms !). The infilled histograms corre- 
spond to the charge states of the fast wind, and the outlined 
histograms correspond to the slow wind. The bottom his- 
togram shows the distribution of mass ratios for the remnant 
grid. Lower Z elements are quickly pushed to high ionization 
states and exhibit a longer recombination time. lonization 
in the slow wind models appears to occur more slowly than 
ionization in the fast wind, but ultimately reaches a higher 
charge state. This is due to the higher density material lead- 
ing to more collisional ionization. 


an emission line's centroid energy and luminosity, as can 
be seen in Figures 14 and 15. Because of this strong den- 
sity dependence, lower mass models can exhibit softer 
centroid energies, especially when the remnant is still 
young. These centroids do still harden in time, and af- 
ter ~ 1000 years they approach the energy which distin- 
guishes SNR Ia from CC SNR, as defined in Yamaguchi 
et al. (2014). 

Examining remnant luminosity as a function of rem- 
nant size yields additional insight into the mass-loss his- 
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Fe-Ka centroids for the fast wind 
(Uwina = 15 kms ^!) models. The centroid luminosities 
are plotted against centroid energy for the: forward shock 
(blue surface), the reverse shock (grey surface), and model 
isochrones for the forward shock (dashed lines). FS values 
correspond to the infilled histograms, and RS values to the 
outlined histograms. 


Figure 14. 


tory of SNR. As can be seen in Figure 16, and Figure 17, 
the Fe-Ka luminosity also correlates well with FS ra- 
dius for a given remnant age. Comparing these results 
to spectral data yield remnants which seem to be rea- 
sonably well constrained by a wind-driven model. As an 
example, the centroid energy, luminosity, and remnant 
size for IC 443 seem to agree with the derived age of 
~ 3000 years and a ZAMS mass of ~ 16 — 25M, from 
Troja et al. (2008). 

Even in cases where remnants do not follow the ex- 
pected behavior of our models, we can still infer some 
information about the mass-loss history. For instance, 
in the case of objects where the centroid energy is softer 
than expected, given its measured size and inferred age, 
we can infer that the Fe in the remnant is not highly 
ionized, which can be interpreted as meaning a lower 
density in the shocked region than is indicated by the 
wind model. Lower than expected Fe-Ka luminosity at 
a given centroid energy can also indicate a lower than 
modeled Fe abundance in the shocked region. 

The current generation of X-ray observatories are ca- 
pable of measuring the Fe-Ka centroid because of its 
luminosity and relative isolation from other lines. Next 
generation observatories, such as XRISM and Athena, 
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Figure 15. Fe-Ka centroids for the slow wind 
(vwind = 5 kms ^!) models. The centroid luminosities 
are plotted against centroid energy for the: forward shock 
(blue surface), the reverse shock (grey surface), and model 
isochrones for the forward shock (dashed lines). FS values 
correspond to the infilled histograms, and RS values to the 
outlined histograms. 


senas 148.3 yrs 411.7 yrs s... 885.0 yrs 15 1922.0 yrs 1 4474.0 yrs 
agus 264.9 yrs 111 607.9 yrs xass 1295.0 yrs 1:1: 2906.0 yrs reire 7000.0 yrs 


Fe-Ka log 1o Luminosity (photons/s) 


0 10 20 30 
Forward Shock Radius (pc) 


Figure 16. Fe-Ka centroid luminosities for the fast wind 
(Uwina = 15 kms !) models. The luminosities are plot- 
ted against remnant size for the total shocked region (blue 
surface) along with model isochrones for the shocked region 
(dashed lines). 
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Figure 17. Fe-Ka centroid luminosities for the slow wind 
(Uwina = 5 kms!) models. The luminosities are plotted 
against remnant size for the total shocked region (blue sur- 
face) along with model isochrones for the shocked region 
(dashed lines). 


will be able to resolve lines in more crowded regions of 
the X-ray spectrum. Many lower-Z elements such as 
Mg or Si still have issues related to either being ionized 
beyond the He-like state, or with other emission lines 
overlaying the transitions, but one possible element that 
is sufficiently separated, while maintaining a significant 
population of He-like ions is S, as discussed in § 4. 


3.2. Hydrogen-like to Helium-like Line Ratios 


Separating the Fe-Ko and Fe-Lya lines is something 
that is within the ability of some current X-ray missions. 
Comparing the relative line strengths of Fe XXV (Fe- 
Ka) and Fe XXVI (Fe-Lya, Fe-Ly8) emission lines can 
be used to probe the ionization state of the remnant. 
In the case of models exhibiting wind-driven mass loss, 
the H- to He-like line ratios also show a dependence on 
the ZAMS to ejecta mass ratio of the model. Higher 
mass ratio models correspond to lower mass-loss rates, 
and lower circumstellar densities. We present the ratios 
of the He- to H-like lines in the shocked circumstellar 
material in Figures 18, and 19. The line ratios in general 
decrease as a function of ZAMS to ejecta mass ratio and 
also as a function of time. 

In Figures 18 and 19, we also include a measured line 
ratio from Koyama et al. (2007) for Sgr-A East. The 
measurement is in the range of our models, although 
it seems to imply a younger than expected remnant. 
The H-like to He-like ratio puts it at ~ 400 — 900 yrs, 
while the Sgr-A East centroid and radius values put it 
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Figure 18. H-like to He-like line ratios for the fast wind 
case. The horizontal dashed line and shaded region corre- 
spond to a measurement of the line ratio for Sgr-A East 
(Koyama et al. 2007). The overall trend is driven by the 
much shorter recombination time and much higher ioniza- 
tion temperature of Fe XXVI compared to Fe XXV. 
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Figure 19. H-like to He-like line ratios for the slow wind 
case. The horizontal dashed line and shaded region corre- 
spond to a measurement of the line ratio for Sgr-A East 
(Koyama et al. 2007). The mechanism and trend are similar 
to the one described in Figure 18, but the increased CSM 
density leads to increased formation of Fe XXVI compared 
to the fast wind case. 


at around ~ 1000 yrs old, which is consistent with the 
literature (Koyama et al. 2007). Sgr-A East seems to be 
over-ionized compared to our models, which also seems 
to agree with the literature. We do note though that it 
is possible Sgr-A East is not a CC SNR, but instead the 
product of a Type Iax SN, which might explain the high 
line ratio (Zhou et al. 2020). 


4. COMPARISONS WITH SNR PROPERTIES 
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While the Fe-K line centroid and luminosity of SNRs 
is a good discriminator between core collapse and Type 
Ia remnants, not all core collapse supernova remnants 
have strong Fe-K line emission. In this section, we ex- 
plore how some of the other SNR spectral and dynami- 
cal properties compare with our models. We focus on a 
handful of relatively young objects which have been pre- 
viously studied in detail with Chandra. We explore how 
the He-like sulfur line centroid can be used to discrim- 
inate between models, as it is typically bright in SNR 
X-ray spectra, even when iron is not. 

We list the handful of objects we explore in Table 1. 
For each SNR, we fit the sulfur He-like triplet to a Gaus- 
sian, the normalization of which is simply the line flux. 
In order to convert the line flux to a luminosity, dis- 
tances to each object are pulled from the literature. 


In total, we examine 2800 model spectra with ages to 
7000 years and radii out to S 30 pc, though we note 
that models which extend to these larger radii may also 
interact with swept-up main sequence wind, which we 
do not consider here. This is a large, but finite set of 
modeled data, and given that simulation produces X-ray 
spectra on a logarithmic timescale, most of our modeled 
spectra are skewed to younger ages. We therefore use 
a Gaussian Mixture Model to represent the distribution 
of our simulated spectra, and fill in the gaps. This is 
discussed in detail below. 

For each modeled spectra, we use XSPEC to simulate a 
Chandra ACIS-S observation of the simulated remnant. 
We then fit the He-like line in the simulated spectrum to 
a Gaussian, in order to get the triplet centroid and lumi- 
nosity. These data are shown in Figure 20, where we plot 
the line centroids and luminosities as a function of for- 
ward shock radius and progenitor ZAMS mass. We also 
overplot the measured spectral characteristics of SNR 
listed in Table 1. 

As seen in Figure 20, the simulations show some clus- 
tering when comparing the spectral properties to one 
another (Fig. 20, left) or when comparing spectral prop- 
erties to dynamical characteristics. Additionally, when 
classifying by ZAMS mass, the lower mass models oc- 
cupy a distinct part of parameter space, while higher 
mass models occupy a separate region. For instance, 
when comparing line centroids to line luminosities, with 
only a few exceptions the models with M > 25 Mo tend 
to lie in the part of parameter space with higher lumi- 
nosities but lower centroids. Likewise, with a handful of 
exceptions, those progenitor models with Mzams S 15 
Mo tend to trend towards lower luminosities, but higher 


line centroids. This trend seems largely independent of 
the remnant age. 

In terms of the line centroids vs radius, clustering in 
the models is a little less clear, but there is a clear spread 
in line centroid at low radius (Fig. 20, right). This clus- 
tering of the model data suggests that the distributions 
of these parameters can be simulated as a mixture of un- 
derlying kernels. We can then plot measured variables 
like those listed in Table 1, without being distracted 
by gaps in the simulated data. To do this, we first as- 
sume that the variables we are comparing, the ZAMS 
mass, model radius, measured sulfur line centroid, and 
measured sulfur line luminosity, can be represented as 
smoothly varying functions. The ZAMS mass is our fun- 
damental input variable, and all other quantities derive 
from it. The assertion that our model grid is drawn 
from a smooth, underlying kernel is reinforced by the 
fact that the SNR ejecta evolution and subsequent ion- 
ization of material are continually varying funtions of 
the SNR age. 

We begin by modeling the distributions of variables 
named above as a mixture of Gaussians. Since we are 
trying to model the underlying distribution which may 
not necessarily be Gaussian, the number we choose can 
be arbitrary. However, we can use the so-called Bayesian 
information criterion (Liddle 2007) to compute the like- 
lihood that the underlying distribution is the sum of a 
number of Gaussians. We iterate on the intial choice 
of underlying density functions, and then compute the 
maximum likelihood that that number of models de- 
scribes the data. In the case of our models, we find 
that the criterion converges at 13 Gaussians. Upon de- 
termining the best number of models to use to describe 
the modeled properties, we replicate the distribution in 
order to see where the data are clustered and how these 
clusters relate to observed SNR properties. 

In Figure 21, we show the same model properties as in 
Figure 20, except now we show the distribution which is 
cloned from our model distribution. As seen in the left 
hand side of Fig. 21, there is definite clustering in the 
model spectral properties. When looking at the cloned 
distribution of forward shock radii vs He-like line cen- 
troid, the clustering becomes more apparent. In par- 
ticular, the higher mass progenitor models tend to have 
lower line centroids than the lower mass models. It is im- 
portant to note that this distribution includes both fast 
and slow winds, and that is revealed in Fig. 21 (right), 
at the smallest radii, where there is a thin strip of high 
mass models at small radius, overlying a broader strip 
of low mass models, while at low centroid energy, there 
is a thin strip of low mass models. Conceptually, this 
suggests that the low mass models which expand into 
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Table 1. Properties of Sulfur He-like triplet for Selected Supernova Remnants 


SNR Distance Res tsnr@ Lx? Centroid Energy Ref. 
kpc pe yr 107? photons s^! keV 

1E 0102-7219 62 6.8 0.2 1738 + 175 17.8 + 3.6 2.445 + 0.009 1,2 
N132D 50 12.5 + 0.2 2450 + 195 173.8 + 15. 2.439 + 0.002 3,4 
Cas A 3.4+0.3 26+ 0.1 340 + 15 7.76 + 0.4 2.449 + 0.001 5 
G11.2-0.3 5.0: 0.5 2.60.3 1400 - 2400 2.88 + 0.5 2.430 + 0.002 6 
Kes 79 71+06 12.4 + 0.8 ~ 6000 0.2 + 0.1 2.435 + 0.007 7 
Kes 73 85+0.8 63+ 0.3 ~ 1800 4.26 + 0.5 2.435 + 0.002 8 
Kes 75 19+2.0 97+ 1.0 ~ 1000 4.30 + 0.75 2.438 + 0.003 9 
RCW 103 3.8 +0.4 50+ 0.5 ~ 4400 0.81 + 0.1 2.431 + 0.002 10 
G292.0+1.8 6 + 0.5 7.7+0.8 3000 — 3400 8.2 + 0.1 2.442 + 0.003 11,12 


"SNR ages are determined either by optical or X-ray proper motion studies of ejecta, or with com- 
parisons to self-similar models. Those with ranges or which are given as approximate values were 


estimated via model comparisons. 


b X-ray luminosity is for the 2.35 — 2.55 keV range. 


References—[1]: Xi et al. (2019); [2]: Banovetz et al. (2021); [3]: Sharda et al. (2020); [4]: Law et al. 
(2020); [5]: Patnaude et al. (2015); [6]: Borkowski et al. (2016); [7]: Sun et al. (2004); [8]: Borkowski 
& Reynolds (2017); [9]: Chevalier (2005); [10]: Braun et al. (2019); [11]: Ghavamian et al. (2005); 


[12]: Bhalerao et al. (2019) 
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Figure 20. Left: He-like sulfur line centroid for our models, colored as a function of progenitor mass. The distribution includes 
both the fast and slow wind cases, and includes models to ages of 7000 years. Right: He-like line centroid vs model forward 
shock radius. In both cases, the data from Table 1 are overplotted. We highlight the measured properties for 1E 0102-7219, 


Kes 73, and G292.0+1.8, which are discussed in the text. 


their low density winds do not efficiently heat and ion- 
ize the ejecta, and conversely the high mass models in 
dense winds quickly ionize the ejecta. 

We conclude this section by quantitatively comparing 
our models to the properties of some of the SNR listed 
in Table 1. We begin with the young O-rich SNR 1E 
0102-7219, located in the Small Magellanic Cloud. By 
measuring the proper motion of several undecelerated 
optical knots, Banovetz et al. (2021) was able to derive 


an age of ~ 1700 yrs, while Xi et al. (2019) estimated a 
progenitor mass of 15-40 Mo, based on the amount of 
circumstellar material they estimate the forward shock 
to have swept up. With a He-like sulfur line luminos- 
ity of 1.8x104 photons s~! and a centroid energy of 
2.445 keV, it lies along the boundary between ZAMS 
models with masses of 15-25 Mo. Comparing its for- 
ward shock radius to its sulfur line centroid (Fig. 21, 
left) also suggests that the progenitor mass was ~ 20 
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Figure 21. Gaussian Mixture Model representation of the simulated supernova remnant properties, colored by ZAMS mass. 
Left: He-like sulfur line centroid versus luminosity. Right: He-like line centroid versus modeled forward shock radius. In both 
cases, the data from Table 1 are overplotted. We highlight the measured properties for 1E 0102-7219, Kes 73, and G292.0--1.8, 


which are discussed in the text. 


Mg. Additionally, the high line centroid energy sug- 
gests that it is interacting with a dense wind. Xi et al. 
(2019) modeled the circumstellar environment inferred 
mass loss rates < 1074 Mo yr-! which, while high, are 
broadly consistent with the rates which come out of the 
red supergiant phase in our models. 

Next, we consider Kes 73. Kes 73 has a similar size 
and age as 1E 0102-7219, but is otherwise quite dif- 
ferent. As seen in Table 1, we measure a low sulfur 
line centroid and luminosity. As seen in Figure 21, its 
spectral and dynamical properties straddle a boundary 
between lower mass progenitor models in a dense wind, 
and higher mass models in a low density wind. Based on 
an analysis of a deep Chandra observation, Borkowski 
& Reynolds (2017) argued for a low mass progenitor 
(Mzams S 20 Mo) which resulted in a "fairly typical” 
IIP supernova. In contrast, Chevalier (2005) compared 
the properties of the remnant to self-similar models, and 
inferred a more massive progenitor which would result 
in a IIL/b SN. Our models suggest a progenitor mass 
of 2 22 Mo, which is in agreement that the remnant is 
from a IIL supernova. 

Finally, we consider the case of G292.0--1.8 (hereafter, 
G292). Located at a distance of = 6 kpc (Gaensler & 
Wallace 2003), G292 has an angular size of ~ 9’, cor- 
responding to a physical diameter of 16 pc (Park et al. 
2007). Based on the proper motion of fast moving opti- 
cal knots, Ghavamian et al. (2005) estimated its age to 
be ~ 3000 yr, consistent with the spin-down age of the 
central pulsar (2900 yr; Camilo et al. 2002). Based on 
the X-ray properties of the shocked circumstellar mate- 
rial, Lee et al. (2010) estimated a mass loss rate of the 
progenitor of M ~ 2-5 x 1075 Mo yr-, for a stellar 
wind speed of 10 km s^ !. 

Recently, Bhalerao et al. (2019) performed an in-depth 
study of a 0.5 Megasecond Chandra ACIS-I observation 


of G292. Based on their detailed spectral analysis and 
estimates for ejecta masses, they suggest that the pro- 
genitor mass of G292 lies within the range of 13 Mọ € 
M S30 Mo. As seen in Table 1, and seen in Figures 20 
(left) and 21 (left), our models suggests a progenitor 
mass of = 20 Mo. Comparisons with the radius vs line 
centroid plots are consistent with this assertion. We 
therefore suggest that G292 is the remnant of a Type IIL 
SN, consistent with previous assert by Chevalier (2005). 

One of the unanswered issues in Patnaude et al. (2015) 
was that while the small grid of models used in that ini- 
tial study could reproduce the observed iron line cen- 
troids and luminosities seen in a large number of rem- 
nants, there was a poor relation in the dynamics between 
observed SNR radius and the spectral properties. Those 
authors asserted that this could possibly be addressed by 
investigating models with larger mass loss rates. They 
argued that the larger mass loss rate would not have 
a large impact on the dynamics, but would impact the 
luminosity and other spectral qualities. As seen in this 
section, as well as in Figures 16 and 17, this hypothesis 
appears reasonable. 


5. APPLICATION TO NEXT GENERATION 
OBSERVATORIES: SPATIALLY UNRESOLVED 
REMNANTS 


Next-generation X-ray observatories such as a XRISM 
and Athena will have unprecedented spectral resolu- 
tion compared to previous instruments. Such resolution 
will allow additional observables such as the S centroids 
mentioned above, or the He-like to H-like ratios of el- 
ements not examined here. While these instruments 
lack the spatial resolution necessary to resolve all but 
the largest and nearest remnants, their integrated spec- 
tra will carry important diagnostic information. For in- 
stance, Doppler shifts and absorption can be used to 
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separate emission from different remnant regions and 
examine the spatial structure of the remnant. We defer 
discussion of these methods to a later publication, and 
choose to focus on what can be determined directly from 
the integrated spectra. 


5.1. SNR Spectral Features for XRISM and Athena 


Next-generation X-ray telescopes will be able to ob- 
serve entire SNRs in the local universe with observa- 
tions of z 10 — 100 kiloseconds. Objects at distances 
of = 10kpc will not be resolvable to instruments such 
as Athena or XRISM, but the integrated spectra still 
contains a great deal of information. The resolving 
power of these new instruments allows for probing the 
He-like and H-like emission from lower-Z elements, as 
can be seen Figure 22. These simulated observations 
show what can be reasonably expected for a middle-aged 
SNR located within the Small Magellanic Cloud (SMC) 
for XRISM, and Athena. He-like transitions are clearly 
resolvable for Mg, Si, and S, as are most H-like lines. 
He-like and H-like Fe emission still exhibit a model de- 


Athena. 
XRISM 


Flux (photons/s/cm?/keV) 


Energy (keV) 


Figure 22. Athena (blue) and XRISM (orange) spectra 
plotted against a 100 eV resolution spectrum (black) for a 
16.7Mọ model in a slow wind normalized to the SMC for a 
300 kilosecond exposure. The vertical lines indicate H-like 
(dashed), and He-like (solid) transitions for Mg (cyan), Si 
(magenta), and S (black). 


pendence, even for the example models presented above. 
Figure 23, presents side-by-side comparisons of XRISM 
and Athena for the three example remnants with a slow 
wind. Even with short exposure times, the He-like lines 
are still easily identified, and the H-like emission is dis- 
cernible in most cases. 


5.2. Helium-like Emission and Line Ratios 


The spectral resolution available to next generation 
observatories will allow for the probing of emission lines 
from individual transitions. Hitomi observations have 
already shown that it is possible to resolve the Forbidden 
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Figure 23. XRISM (left), and Athena spectra plotted 
against the original model spectra for the total shocked re- 
gion of the three example models, normalized to the SMC. 
The XRISM spectra were exposed for 300ks, and the Athena 
spectra were limited by the large number of counts at ~ 
1 keV and could only be exposed for 3ks, 20ks and 300ks 
(top to bottom). The He-like and H-like lines of Fe are indi- 
vidually resolvable with both instruments. 


(F,z), Intercombination (I, x, and Ig, y), and Resonance 
(R,w), transitions of He-like Fe (Hitomi Collaboration 
et al. 2016). The ratios of these lines can be used to 
determine the underlying physical conditions present in 
the plasma (Porquet et al. 2010). The g-ratio is defined 


as pet 
aja to (2) 


w 
and has a dependence on the electron temperature. Ad- 
ditionally, the r-ratio is defined as 

r-c-y 


r (ne) = (3) 


z 


and has a dependence on the electron number density. 
Unfortunately, there is not sufficient difference in the 
dependence on thermodynamic variables to differentiate 
between models, as they all represent similar NEI con- 
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Figure 24. The Fe G-ratio extracted from the integrated spectra for the shocked region. In general, a larger ZAMS mass 


results in a lower intial G-ratio, with values eventually converging to ~ 4.25 (horizontal line). 


Lower mass models exhibit an 


early time peak, along with those models that experienced an early time blue loop (15.7Mo, 16.1Mọ) 


ditions. Additionally, the r-ratio for high-Z elements 
requires significantly larger densities than are present in 
these SNR models. 

'The g-ratio does exhibit time and mass-dependent be- 
havior, with the ratios of all models converging to + 4 
within ~ 10? yrs. Additionally, models that exhibited a 
low mass loss rate (high mass ratio) often started with 
a higher initial g, and exhibited a peak within the first 
~ 10? yrs, before settling near the limit. Objects where 
the mass loss was much lower than would usually be 
expected (ie. the blue loops) exhibited a lower initial 
g-ratio, which peaks within ~ 100 yrs before quickly re- 
turning to the limit. Meanwhile, models that exhibited 
a higher mass loss rate started with a significantly lower 
g-ratio, and rose to the limit. Figure 24 explores these 
behaviors in detail. There is a great deal of variability, 
especially at very early times when the total fraction of 
shocked material is still low, but when considering rem- 
nants younger than ~ 100 — 200 yrs, He-like line ratios 
could be a useful probe of the mass loss rate, and the 
ZAMS mass of the progenitor. As in the forward shock 
only case, the line ratio begins by decreasing, reaching 
some minimum value before rising again. 

'The H-like to He-like line ratio behavior also changes 
slightly when considering emission from the entire 
shocked region. Unlike in the forward shock only case, 
the line ratio begins by decreasing, reaching some min- 
imum value before rising again. This is consistent with 
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Figure 25. The ay line ratios for the integrated shock 


region as a function of time and mass ratio for the slow 
wind (top) and fast wind (bottom). The early time behav- 
ior is largely the same, but there is a second rise at times 
(2 150 yrs). 


what is described in Patnaude et al. (2015), where He- 
and H-like line emission is first dominated by shocked 
circumstellar material. But as the SNR evolves, shocked 
ejecta becomes more important to the integrated X-ray 
emission, causing it to rise, and altering the line ratios 
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in ways that are dependent upon both the ejecta mass 
and circumstellar environment. Overall, the trend of 
line ratio vs. mass ratio still exists as in the FS case, 
with lower mass ratio models exhibiting a higher line 
ratio for a given remnant age. There also exists a flat- 
tening of the data below a mass ratio of 0.7 in the fast 
wind, with those models maintaining a minimum line 
ratio of z 0.03 before rising again. The division appears 
more strongly in the fast wind case, and appears to be 
washed out by the increased density of the slow wind. 
These trends can be seen in Figure 25. 


6. SUMMARY AND FUTURE WORK 


Supernova remnant evolution is dependent on the 
mass loss history of the progenitor object. When consid- 
ering line-driven winds it is possible to discern between 
progenitor models based on the spectra of the resulting 
remnant. Ejecta masses tend to cluster around a narrow 
range, with the main difference in composition being the 
CO core. There is variation in mass loss based on the 
wind prescription used, which leads to variation in the 
CSM density, but the overall trends are similar within a 
given wind prescription, and show an identifiable trend 
based on progenitor mass. the Fe-Ka centroid is a rea- 
sonable discriminant of the progenitor type (Ia or CC) 
and the ZAMS mass. The He-like to H-like line ratios 
are a good proxy for the ionization state of the remnant. 

In cases where the stellar mass loss is well described 
by an isotropic wind, remnant information can be re- 
liably recovered using the Fe-Ka centroid and the ra- 
dius. Without examining the individual composition 
and metallicity of a remnant, it is still possible to exam- 
ine some aspects of the X-ray spectra and make infer- 
ences about the mass range, as in the case of G292.0+1.8 
and others. Future X-ray missions will allow for addi- 
tional observables such as the He-like line ratios which 
can be used to probe the CSM density and progenitor 
mass even in the case of a spatially unresolved remnant. 

Models that have a more complex mass loss history, 
such as the blue-loop models will lead to a more com- 
plex CSM, as the lower mass loss rate during the loop 
will lead to a low density shell embedded within the 
denser RSG wind. That aspect is currently neglected, 
along with the high speed winds expected from YSG, 
BSG, and Wolf-Rayet stars. A more detailed investiga- 


tion of those remnants is warranted. The current mass 
loss prescriptions also led to partially stripped envelope 
SNe only for higher mass models (M 2 22M). This 
is broadly consistent with the progression of mass strip- 
ping in supernova progenitors when considering IIP vs 
IIL SNe. For lower mass stars, additional mass loss pre- 
scriptions are required to remove the hydrogen envelope 
such as wave-driven mass loss or binary interactions. 

Finally, Projection effects combined with shock 
Doppler boosting and intra-remnant absorption can be 
used to extract additional information from spatially un- 
resolved remnant spectra. Modeling of these effects is 
currently being incorporated into the above framework, 
and will be the subject of a future work. 
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